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We discuss the model OaniQ^) recently proposed for the QCD running coupling as{Q^) in the Euclidean domain 
on the basis of the "asymptotic-freedom" expression and on causality condition in the form of the Q^-analyticity. 

The model contains no adjustable parameters and obeys the important features: (i) Finite ghost-free behavior 
in the "low Q^" region and correspondence with the standard RG-summed UV expressions; (ii) The universal 
hmiting value Qan(O) expressed only via group symmetry factors. This value as well as the San behavior in the 
whole IR region < turns out to be stable with respect to higher loop corrections; (iii) Coherence between 
observed as{M^) value and integral information on the IR Qs(Q^) behavior extracted from jet physics. 



1. INTRODUCTION 

This presentation is devoted to the review and 
discussion of a new analytized model expression 
aan(Q^) for the QCD running coupling recently 
devised Q by combining the RG-summed expres- 
sion ctsiQ^) with the demand of analyticity, that 
is causality, in the complex plane. This proce- 
dure "cures" the IR ghost-pole trouble by an ad- 
ditional contribution that is non-analytic in the 
coupling constant and at the same time preserves 
the asymptotic freedom property and correspon- 
dence with perturbation theory in the UV . 

The "analytization procedure" elaborated in 
the mid-fifties (see Ref.g) consists of three steps: 

- Find an explicit expression for ajidQ^) in 
the Euclidean region > by standard RG 
improvement of a perturbative input. 

- Perform the straightforward analytical con- 
tinuation of this expression into the Minkowskian 
region ReQ'^ < 0,lmQ^ = — e. Calculate its 
imaginary part and define the spectral density by 
PRoia; a) ^ lma_RG(-o- - ie, a). 

- Using the spectral representation of the 
Kallen-Lehmann type with p^c in the integrand, 
define an "analytically-improved" running cou- 
pling aan(Q^) in the Euclidean region. 

Being applied to a(Q^) in the one- and two- 
loop QED {or QCD} case, this procedure pro- 
duces (see Ref. |] {or [|l]||}) an expression 
cian(Q^) with the properties: 

(a) it has no ghost pole, 

(b) in the complex a-plane at the point a = it 



possesses an essential singularity ^ exp(— l/a/3i), 
with the one- loop beta-function coefficient, 

(c) in the vicinity of this singularity for real 
and positive a it admits a power expansion that 
coincides with the perturbation input, 

(d) it has the finite UV {or IR} hmit Att/Pi 
that does not depend on the experimental value 
a ~ 1/137 {or A}. 

In the one-loop QCD case this expression is of 
the form 

^"^^ ' lnQ2/A2 + A2-Q2 

with a = /3ias/47r and /3i = 11 — (2/3)n/. 

The same procedure being applied to the two- 
loop case yields |||,|l|,^ a more complicated expres- 
sion with the same essential features (see, e.g., 
Eqs. (I) - d) below). 

2. THE MODEL FOR QCD RUNNING 
COUPLING 

The "analytic" coupling constant, Eq. (|l|), in- 
stead of ghost pole has a weak singularity at 
Q2 ^ and its IR limit a^^O) = An /Pi depends 
only on group factors. Numerically, for n/ = 3, we 

have ain^(O) = 47r/9 ~ 1.4. 

Note, to relate A, the QCD scale parameter, 
to a(/i^), in our case we have to change its usual 
one-loop definition for A^ = n"^ exp[— (/)(a(/i^))] 
with the function (j) defined by the transcedental 
relation 



a=lMa)-l/[exp(0(a))-l] 



(2) 
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which is consistent with the symmetry property 

c^{a) = -0(1 - a) . (3) 

The corresponding beta-function 

M \ 1 I exp0(a) 

P(a) = ~T^TT^ H 9 (4) 

02(a) [exp0(a)-l]2 

/3(1 — a) and obeys the 



with 



is also symmetric /3(a) 
second-order zero at a 



1 - see Fig. 1. 
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Figure 1. One- loop /3(a) function. 



Usually, in perturbative QFT practice, we are 
accustomed to the idea that theory supplies us 
with a set of possible curves for as{Q) and one 
has to fix the "physical one" by comparing with 
experiment. Here, Eq. (^) describes a family of 
such curves for aan(Q2) forming a bundle with 
the same common limit at = 0. 

For the two-loop case, we start with the run- 
ning coupling written down in the form 
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; + 51n(l + //6) 



I ^ In 



A2 



(5) 



This 



with h = (32/{Piy and (32 = 102 - 38/3n/. 
expression corresponds to the result of exact in- 
tegration of the two-loop differential RG equa- 
tion explicitly resolved by iteration. It generates 
a popular two-loop formula with the Inl/P term. 
For the spectral density, we have 



I(L) 



R^L) + P{L) ' 



L = In 



A2 



(6) 



R{L) ^L + b\n ^{l + L/b)^ + {n/by , 
1 (L) = n + b arccos 



^J{b + Lf+TT^ 



(7) 
(8) 



Now, to obtain a'^^iQ'^), 



one has to substitute 



Eq. 



into the r.h.s. of 



da- 



+ Q2 - 



(9) 



The one- loop result, Eq.(|I|), follows from Eqs. (g) 
- (H) at 6 = 0. However, in the two-loop case, the 
integral expression thus obtained is too compli- 
cated for being presented in an explicit form. For 
a quantitative discussion one has to use numerical 
calculation. 

Nevertheless, for a particular value at = g 
we can make two important statements. First, 
the IR limiting coupling value Q;an(0) does not de- 
pend on the scale parameter A. Second, this value 
turns out to be defined by the one-loop approxi- 
mation, i.e., its two-loop value coincides with the 
one- loop one (for detail see ref. Q). 

Thus, the Q;an(0) value, due to the RG invari- 
ance, is independent of A and, due to the analytic 
properties, of higher corrections. This means that 
the analyticity stabilizes the running coupling be- 
havior in the IR, makes it universal. 

Moreover, the whole shape of the aa.n{Q^) evo- 
lution turns out to be quite stable with respect 
to higher corrections. The point is that the uni- 
versality of aan(O) gives rise to stability of the 
oiniQ^) behavior with respect to higher correc- 
tion in the whole IR region (at ~ A^). On the 
other hand, stability in the UV domain is reflec- 
tion of the property of asymptotic freedom. As 
a result, our analytical model obeys approximate 
"higher loops stability" in the whole Euclidean 
region. Numerical calculation (performed in the 
MS scheme for £ = 1, 2 and 3 cases with n/ = 3) 
reveals that ain (Q^) differs from ain {Q^) by less 
than 10% and ain (Q^) from ain{Q'^) within the 
1% limit, as it is shown in Fig. 2. 

The idea that QCD running coupling can be 
frozen or finite at small momenta was considered 
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Figure 2. "Higher loop stability" of the analytic 
solution. The normalization point is on the t 
lepton scale: as{M^) = 0.34. 



in many papers (see, e.g., the discussion in Ref. Q 
and interesting theoretical model in Refs. |^|). 
There seems to be experimental evidence in fa- 
vor of such behavior. As an appropriate object, 
for comparison with our construction, we use the 
average 
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AiQ) ,k 



(10) 



that people manage to extract from jet physics 
data. Empirically, it has been claimed 1^, 
that this integral at Q — 1-^2 GeV turns out 
to be a fit-invariant quantity with the estimate: 
A{2 GeV) = 0.52 ± 0.10. Numerical results for it 

-(2) . A 

obtained by the substitution aan into Eq. (10) 
with Aan — 710 ±110 MeV, corresponding to 
Q;an(M^) = 0.36 ± 0.02, are summarized in the 
Table. 



aa„(M^) 


0.34 


0.36 


0.38 


Ai'J(MeV) 


610 


710 


820 


A(2)(2GeV) 


0.48 


0.50 


0.52 



Note that a nonperturbative contribution, like 
the second term in the r.h.s. of Eq. (|l|), reveals it- 
self even at moderate Q values by "slowing down" 
the velocity of as{Q) evolution. For instance, 
at (5 = 3 GeV it contributes about 4%, which 
could be essential for the resolution of the "dis- 
crepancy" between "low-Q" and direct Zq data 
for asiMz). 



As far as we have no explicit expression for the 
Q!an('9^) in the two-loop case, for a qualitative dis- 
cussion we can use an approximate formula pro- 
posed in Ref. |^ which can be written in the form 
of Eq. (m with the substitution 
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A2 



exp 



In ^ + In \ / In^ ^ 
A2 V A2 



-47r2 



(11) 



With appropriate redefinition of A, the accu- 
racy of this approximation, for Q > A, is no less 
than 5%. At the same time, it produces only a 
3% error in the A{2 GeV) value. 

3. DISCUSSION 

It is important to discuss the possibility 
of using aan(Q^) in multi-prong QFT objects 
Fffc) (si, Sfc) and, particularly, in observables 
M{...Qi, ...,qj...) with some arguments — 
—Qj > being time-like (Minkowskian) and 
some others fixed on the mass shell. 

Here, we have in mind a few different items: 



1 . The possibility of using RG for fc-prong ver- 
tices F(fc)(si, Sfc) , fc > 3; 

2. The technology of using Q:an(Q^), originally 
defined for Q| > 0, in observables 
M{...Qj, (?j ...) with time-like arguments; 

3. The expediency of the aan(Q^) continuation 
into the Minkowskian region < 0; 

4. The scheme dependence of analytic running 
coupling. Relation to divergencies. 

Our preliminary comments on these issues are: 
F As it is well known from the old investiga- 
tions, the use of RG, rigorously deduced from 
Dyson renormalization transformations (with fi- 
nite real counterterm coefficients Zi), is justified 
only in the Euclidean domain and involves a si- 
multaneous change of a scale for all arguments of 
F(j.)(sj), i.e., Sj s'j — tsj. This restricts the 
possibility of UV analysis by the so-called non- 
exceptional momenta. 

H. Nethertheless, in some special cases it 
turns out to be possible to apply the RG tech- 
nique to analyse multi-prong object by involving 
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some additional means like spectral representa- 
tions light-cone expansion or amplitude fac- 
torization 

In any of above-mentioned cases the RG proce- 
dure results in solving group equations for real 
functions, like spectral densities and light-cone 
expansion coefficients. This solution involves only 
real running coupling aan(Q^) values for Eu- 
clidean arguments. Any observable A4 with time- 
like value of kinematic invariant should be 
treated separately with an adequate procedure of 
analytic continuation of the observable under con- 
sideration. 

III. Due to the last reason, the analytic con- 
tinuation of the running coupling itself, e.g., dis- 
cussing of ttanCQ^) properties in the Minkowskian 
region < 0, in our opinion, has no direct phys- 
ical meaning. 

IV. The last but not least is the property of 
the scheme dependence of the model discussed. 
Formally, in our final expressions, Eqs. 

there is no room for such dependence. This is 
related to the absence of UV infinities with their 
subtraction and renormalization ambiguities. 

The analytical model has an important prop- 
erty. By construction it is free from UV divergen- 
cies. The log squared in the spectral function 
denominator (|^) provides us with convergence of 
non-subtracted spectral representation (||). At 
the very end, it contains only one parameter A 
that has to be defined from experiment. How- 
ever, this needs an adequate procedure (men- 
tioned above in the comment II) of analytiza- 
tion with possible liberating of infinities by con- 
tribution non-analytic in a for the observable con- 
fronted with data. 

To summarize, it can be said that to get more 
satisfactory answers and, correspondingly, more 
complete understanding, it is necessary to con- 
tinue investigation of all four issues. 



edged. 

REFERENCES 

1 D.V. Shirkov and I.L. So lovtsov, JINR R apid 
Comm., No. 2[76]-96, 5, lhep-ph/960436^ . 

2 N.N. Bogohubov, A. A. Logunov and D.V. 
Shirkov, Sov. Phys. JETP 37(10) (1960) 574. 

3 D.V. Shirkov and I.L. Solovtsov, hep- 



ph/9704363| ; Phys. Rev. Lett. 79 (1997) 1209. 

4 A.C. Mattingly and P.M. Stevenson, Phys. 
Rev. D49 (1994) 437. 

5 Yu.A. Simonov, Yad. Fizika, 58 (1995) 1139; 
A.M. Badalian and Yu.A. Simonov, ibid. 60 
(1997) 714. 

6 Yu.L. Dokshitzer and B.R. Webber, Phys. 
Lett. B 352 (1995) 451. 

7 Yu.L. Dokshitzer, V.A.Khoze and S.I. Troyan, 
Phys. Rev. D53 (1996) 89. 

8 I.F. Ginzburg and D.V. Shirkov, Sov. Phys. 
JETP 22 (1966) 234. 

9 A.V. Efrcmov, A.V., Radiushkin, Riv. Nuovo 
Cini. 3 (1980) 1-87. 

Discussions 
K.Chetyrkin 

What do you think about an experimental test- 
ing of different predictions for the IR behavior of 

as{Q) ? 

D. Shirkov 

As I showed in the Table, our model nicely 
correlates the as(MT-) measured value with the 
Khoze- Dokshitzer integral estimate extracted from 
the jet physics. This can be compared, e.g., 
with the A{2 GeV) value corresponding to the 
Badalian- Simonov model |^ which is of an order 
of 0.35. 
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